Frequency Scales for Current Statistics of Mesoscopic Conductors 
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We calculate the third cumulant of current in a chaotic cavity with contacts of arbitrary trans- 
parency as a function of frequency. Its frequency dependence drastically differs from that of the 
conventional noise. In addition to a dispersion at the inverse RC time characteristic of charge re- 
laxation, it has a low-frequency dispersion at the inverse dwell time of electrons in the cavity. This 
effect is suppressed if both contacts have either large or small transparencies. 
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What is the characteristic time scale of the dynamics of 
electrical transport in normal-metal mesoscopic systems 
if their size is larger than the screening length? The most 
reasonable answer is that this dynamics is governed by 
the RC time of the system, which describes the relaxation 
of piled up charge. It is this time that characterizes the 
admittance of a mesoscopic capacitor, 1 the impedance of 
a diffusive contact, 2 current noise in diffusive contacts, 3 
and charge noise in quantum point contacts and chaotic 
cavities. 4 ' 5 The free motion of electrons is superseded by 
the effects of charge screening and therefore the dwell 
time of a free electron is not seen in these quantities. 
The only known exception is the frequency dispersion of 
the weak-localization correction to the conductance. 6 ' 7 
Because typical charge-relaxation times are very short 
for good conductors, the experimental observation of the 
dispersion of current noises is difficult. 8 

Very recently, the third cumulant of current has 
been measured for tunnel junctions. 9 In this case, the 
frequency dispersion is due only to the measurement 
circuit. 10 In contrast to tunnel junctions, chaotic cavi- 
ties have internal dynamics. It is the purpose of this 
paper to show that mesoscopic systems with such dy- 
namics may exhibit an additional low-frequency disper- 
sion in the third cumulant of current, which is absent 
for the first and second cumulants. Physically, this dis- 
persion is due to slow, charge-neutral fluctuations of the 
distribution function in the interior of the system. Such 
charge-neutral fluctuations are similar to fluctuations of 
the effective electron temperature at a constant chemical 
potential. They do not contribute directly to the mea- 
surable current, but they modulate the intensity of noise 
and hence contribute to higher cumulants of current. 11 

The zero-frequency electric noise in open chaotic 
cavities has been calculated by Jalabert et al. 12 us- 
ing random-matrix theory. More recently, these ex- 
pressions were derived semiclassically by Blanter and 
Sukhorukov. 13 Shot-noise measurements on chaotic cavi- 
ties were performed by Oberholzer et al. 14 The third and 
fourth cumulants of current were obtained by Blanter et 
al. 15 for open cavities and in Ref. 16 for cavities with con- 
tacts of arbitrary transparency. Recently, the statistics 
of charge in a cavity was analyzed by two of us. 17 




FIG. 1. Left panel: A chaotic cavity is biased from two 
leads and capacitively coupled to a gate. Right panel: 
Charged fluctuations of the distribution function / decay at 
the short RC time and contribute to the current noise. Un- 
charged fluctuations decay at a much longer time to and do 
not contribute directly to the current noise but result in a 
low- frequency dispersion of the third cumulant of current. 



The system we investigate is shown in Fig. 1. A 
chaotic cavity consists of a metallic mesoscopic conduc- 
tor of irregular shape connected to leads L, R via quan- 
tum point contacts with conductances Gl,r. and trans- 
parencies Tl,r- We assume that G l ,r ^> e 2 /h and the 
bias is much larger than the temperature and frequen- 
cies. The dwell time of electrons in the cavity is given 
by t d — 6 2 RqN f , where Rq = (Gl + Gr)^ 1 is the 
charge-relaxation resistance 7 and Np is the density of 
states in the cavity, which is assumed to be continuous 
and constant. We assume that the dwell time is short as 
compared to inelastic scattering times. Then the mean 
occupation function may be written as a weighted aver- 
age / (e) = RQ{G L fh{() + G R f R {e)) of the Fermi oc- 
cupation factors /l,r of the leads which are step func- 
tions at zero temperature. In the absence of electrostatic 
interaction, fluctuations around /o would relax on the 
time scale Tr> However, due to screening there is a sec- 
ond shorter time scale tq, which describes fluctuations 
of the charge in the cavity. It is given by tq — RqC^, 
where C^ 1 = C _1 + (e^Np)^ 1 is the electrochemical 
capacitance 1 and C is the geometric capacitance. 

Our semiclassical calculations are based on a large sep- 
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aration between the time scales describing the fast fluc- 
tuations of current in isolated contacts and slow fluctu- 
ations of the electron distribution in the cavity. 16 This 
allows us to consider the contacts as independent gen- 
erators of white noise, whose intensity is determined 
by the instantaneous distribution function of electrons 
in the cavity. Based on this time-scale separation, a 
recursive expansion of higher cumulants of semiclassi- 
cal quantities in terms of their lower-order cumulants 
was developed. 11 ' 16 Very recently, such recursive relations 
were obtained as a saddle-point expansion of a stochastic 
path integral. 18 Here we follow the approach of Ref. 18 

Our aim is to obtain the statistics of time-dependent 
fluctuations of the current /z,(t) flowing through the 
left contact of the cavity. In general, the probability 
V[A{t)\DA of a time-dependent stochastic variable A(t) 
(t G [0,T]) is described by the characteristic functional 
S[xa] evaluated with an imaginary field 



V[A(t)] = f VxAe^fo ' dtA XA+ s[ XA ~i XA ] 



(1) 



Functional derivatives S n /(SxA(ti) ■ ■ ■ SxA{t n )) of the 
characteristic functional S[xa] yield irreducible correla- 
tion functions (A(ti) . . . A(t n )). 

To obtain the characteristic functional S[xl] describ- 
ing the fluctuations of we proceed in two steps. 
First, we consider the point contacts (i — L,R) as 
sources of white noise that depend on two common time- 
dependent parameters, the electron occupation function 
of the cavity f(e — eU) and the electrostatic potential of 
the cavity U. Their characteristic functionals are given 
by Si[xi,e] = JdtJ deHi(xi,e) with 
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The correlators (/™} e must be taken from a quantum- 
mechanical calculation 19 ' 20 



ln{l + Tifi(e)[l - f(e - eU)]{e~^ - 1) 
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+rj(e-eU)[l-Me)}(e^~l)}\ x= o. 



(3) 



In a second step, we take into account that the occupa- 
tion function f(e — eU) and the potential U are not free 
parameters but are fixed by the kinetic equation 
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We introduced the fluctuating part of the occupation 
function Sf(e) = f(e — eU) — ,fo(e). The two conserva- 
tion laws are expressed by path integrals over Lagrange 



multipliers A e ,£ and integrated over the fluctuations of 
occupation function and potential to obtain the following 
result for the generating functional 
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x exp{S L [i\ e % X l\ + S R [i\ £ + i$ + S c } , (5) 

where the conservation laws are expressed by the follow- 
ing dynamical action 

S c = -ij^ dtfe CU + J dej5f(e) 



, i d 1 \ P . ,r dU df 



-Glxl J de[f(e-eU)-f L (e)]}. 



(6) 



In the semiclassical regime, this path integral may be 
evaluated in the saddle-point approximation. 18 The sad- 
dle point equations arc nonlinear differential equations 
for the four fields f(e — eU), SU, X e , £ that contain the ex- 
ternal statistical field xl in inhomogeneous source terms. 
They describe the non-linear response of internal fields to 
these sources. In this publication, we are interested in the 
second and third order correlation functions of the cur- 
rent II- For this purpose it is sufficient to expand the 
sum i] = A c + £ = ?7i + 772 /2 + . . . up to second order 
in the external field (r] n oc x£). The result can then be 
substituted into the action (5) to obtain the frequency 
dispersion of the correlation functions. 

We first briefly discuss the frequency dependence of 
current noise to connect our theory to earlier results. 4 
To this end we collect the second order contributions of 
Eq. (5) and eliminate time derivatives using the saddle 
point equations 



S2[X1 



dtde 



(7) 



The zero subscript indicates that the correlators (I") e are 
evaluated at /(e — eU) = /o(e)- We choose a large time 
interval [0, T] and neglect any transient effects. It is then 
most convenient to work completely in Fourier space. 
The linear response of the internal fields (1 + iujTQ)r\\ = 
GlRqXl is entirely governed by the RC time, therefore 
the correlation function (Il{uj)Il(lo')) = 2-k5{uj + w')P 2 
takes a simple form 



P 2 = \Z Q {io)\ 2 / de 



(7E)o|GkM| 2 + (/ 2 ) G| 



(8) 



where we introduced Zq(u>) = (Gl + Gr + iojC^) 1 and 
Gr(lo) — Gr + iivC^. At low frequencies, the current 
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FIG. 2. Cascade expansion for the third cumulant of cur- 
rent, (a) Second cumulant, (b) minimal-correlation third cu- 
mulant, and (c) cascade correction to the third cumulant. The 
external ends correspond to current fluctuations at different 
moments of time and the dashed line, to fluctuations of the 
distribution function. The full circle and triangle correspond 
to the second and third cumulants of white noise generated at 
the contacts. The white triangle corresponds to the functional 
derivative of the second cumulant. 
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FIG. 3. A 3D plot of ReP 3 (wi,w 2 ) for Gl/Gr = 1/2, 
Tl = Tj} = 3/4, tq = 1/3, and to = 10 (dimensionless 
units). 



noise shows correlations between left and right point con- 
tact. At high frequencies, these correlations disappear 
and we observe the bare noise of the left point contact. 
The transition frequency is given by the RC time. To 
complete this result, it remains to evaluate the bare noise 
correlators using Eq. (3). 

We now turn to the much more complicated frequency 
dependence of the third order correlator. For the action 
(6) we find 



dtde 



{{llUrH-XL? + {P R )or,l 
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(iDoim - xl) + (i 2 R )om] vi} 



(9) 



This correlator contains two contributions presented by 
diagrams in Fig. 2. The first one represented by diagram 
b is the minimal correlation result 13 that depends only 
on the RC time. The second is represented by diagram c 
and gives the cascade correction, 11 ' 16 which contains the 
low-frequency dispersion. The equations of motion are of 
the form 

(1 + iujT D )X 2 = R Q A L ( m - xl) 2 + RqArtjI - 
(1 + ilut q )& = R Q B L ( m - xl) 2 + RqBrvI (10) 

and depend on both time constants tq and tjj that 
describe the decay of charged and charge neutral fluc- 
tuations of the occupation function /(e — ell). In 
turn, these fluctuations act back on the current noise. 
We used here a notation Ai = d(lf)o/df, and Bi = 
(Cp/C) J de(—dfo/de)Ai. It is now straightforward to 
substitute the second order response (10) into the ac- 
tion (9). Functional derivatives with respect to the ex- 
ternal field xl then yield the irreducible part of the 
third order correlation function (Il{^i)Il{^2)Il{^3)) = 
2ir6(uji + u>2 + lo 3 )P 3 . The same results are obtained by 
using a recursive diagrammatic expansion of the third 
cumulant. 11 ' 16 



In what follows, we will be interested in the most typi- 
cal case where the dwell time td is much larger than the 
charge relaxation time tq = C'^Rq. Even for this case 
the general expression for P 3 is too long and we present 
here only its limiting values. In the low- frequency limit 
uo t < (C^Rq)- 1 , P 3 is of the form 21 



P 3 (wi,w 2 ) =e 2 l{3G L G R 
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which in general suggests a strong dispersion of the third 
cumulant of noise at oji_ 2 ~ 1/rjj. This dispersion van- 
ishes for symmetric cavities and cavities with two tunnel 
or two ballistic contacts. For a cavity with two tunnel 
contacts, the white-noise sources (3) are linear function- 
als of the distribution function and hence are not affected 
by charge-neutral fluctuations. For the case of two bal- 
listic contacts, Eq. (2) depends only on / and not on 
/j. Then the low- frequency dispersion does not show 
up either because fluctuations of current and distribu- 
tion function are uncorrelated due to the symmetry of 
Hi. Note that Eq. (11) is symmetric with respect to in- 
dices L and R because pile-up of charge in the cavity is 
forbidden at low frequencies. 

In general, P$(u)\,u)2) exhibits a complicated behavior 
(see Fig. 3). It has also a dispersion at the inverse RC 
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FIG. 4. P 3 (w,0) as a function of w for Gl/Ga = 1, 
= 1, Tfj = 0, and to = 10 (dimensionless units). The 
solid, dashed, and dash-dot curves correspond to tq = 1/2, 
tq = 3, and to the case of weak electrostatic coupling tq = td 
(G = co). 

time and additional peculiarities at the scale r^ 1 if one 
of the frequencies or their sum tends to zero. The shape 
of P 3 (cj 1 ,w 2 ) essentially depends on the parameters of 
the contacts. In particular, for a cavity with one tun- 
nel and one ballistic contact with equal conductances 
Gl = Gr = G it exhibits a non-monotonic behavior 
as one goes from oj\ = ll>2 = to high frequencies. A rel- 
atively simple analytical expression for this case may be 
obtained if Try » tq and one of the frequencies is zero: 



1 „ 1 + 2r z n uj' z + T^r^ 4 
P a( „,0) = --^/ (i + ^ )(i+ °«, )2 . ,12, 

The Ps(u), 0) curve shows a clear minimum at lo ~ 
(t>tq) _1 / 2 and the amplitude of its variation tends to 
P 3 (0,0) as t q /t d -» (see Fig. 4). 

The fluctuations of the current in the left-hand con- 
tact of the cavity can be measured as fluctuations of 
voltage across a small resistor attached to it. Based 
on the parameters of a chaotic cavity used in shot-noise 
experiments, 1 our estimates give an inverse dwell time 
of the order of 10 GHz, which is well within the exper- 
imental range for measuring the frequency dependence 
of noise. 8 We also believe that similar low-frequency dis- 
persion of the third cumulant may be observed in other 
scmiclassical systems like diffusive wires, where the dwell 
time is of the same order or larger. Hence our results are 
of direct experimental interest. 

In summary, we have shown that the third cumulant 
of current in mesoscopic systems may exhibit a strong 
dispersion at frequencies much smaller than the charge- 
relaxation time of the system. This effect has a purely 



classical origin and the variations of the cumulant may 
be of the order of its zero-frequency value even if the 
number of quantum channels in the system is large. 
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